日 干 || 


ChinaXiv 合 作 期 刊 


一 种 基于 相位 条 纹 的 高 精度 码 相位 测量 方法 
和 孙 文 杰 "”， 王 兆 瑞 

(1. 中 国 科 学 院 国 家 天 文 台 , 北京 100101; 2. 中 国 科学 院 大 学 ， 天 文 与 空间 科学 学 院 , 北京 
100049 ) 

摘 要 : GNSS (Global Navigation Satellite System， 全 球 导 航 卫 星系 统 ) 接收 机 在 接 
收 到 卫星 信号 后 ， 经 过 捕获 、 跟 踪 ， 得 到 码 相位 与 载波 相位 两 个 基本 的 测量 值 。 然 后 通过 测 
量 值 计 算 卫 星 与 接收 机 之 间 的 距离 。 载 波 相位 测量 值 比 码 相 位 测量 值 具有 更 高 的 测量 精度 ， 
0 量 值 具有 计算 整 周 模糊 度 的 问题 ,在 实际 的 计算 当中 得 到 高 精度 的 相位 测 

结果 需要 付出 更 多 的 计算 代价 。 因 此 ,在 实际 应 用 中 ,得 到 更 高 精度 的 码 相位 测量 结果 有 具 
A 本 文 在 基于 相位 条 纹 的 技术 基础 上 ， 通 过 提取 互相 关 功 率 谱 ， 提 出 一 种 更 
高 精度 的 码 相 位 测量 方法 , 得 到 更 高 精度 的 码 相 位 测量 值 ， 从 而 得 到 更 高 精度 的 码 伪 距 | 测量 
结果 ， 并 通过 实际 的 GNSS 的 信号 验证 了 本 方法 的 有 效 性 
关键 词 : GNSS; 相位 条 纹 法 ; 码 相位 
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卫星 导航 接收 机 通常 使 用 载波 辅助 的 载波 跟踪 环 和 码 跟 踪 环 实现 载波 跟踪 与 码 跟踪 5 。 


其 中 码 跟踪 环 通常 使 用 非 相 干 延 迟 锁定 环 DLL (Delayed Lock Loop) 来 实现 ， 通 过 接收 到 


的 I、Q 两 路 信号 ,分别 产 生 超 前 、 即 时 、 灌 后 的 6 组 信号 ， 利用 GNSS 伪 随 机 码 的 自 相关 特 


性 ， 用 鉴 相 器 计算 出 本 地 产生 的 码 与 接收 到 的 码 的 相位 差 ; 在 获得 相位 差 的 同时 ,将 鉴 相 结 


果 输 入 给 数字 控制 振荡 器 NCO (Numerically Controlled 0scillator)， 调 节 本 地 产生 的 伪 


码 的 相位 ， 使 之 与 输入 的 码 的 相位 进行 对 齐 ， 完 成 码 跟 踪 。 同 时 也 得 到 码 相位 测量 值 ， 可 以 


用 于 计算 伪 距 。 以 GPS (Global Positioning System) 为 例 ，GPS 的 码 速率 为 1. 023MHz， 


一 个 码 片 的 长 度 为 293m， 码 伪 距 测量 误差 约 为 一 个 码 片 长 度 的 1%， 大 约 为 3m 的 误差 。 


为 得 到 更 高 精度 的 码 相位 观测 值 , 国内 外 许多 学 者 尝试 改善 码 环 的 结构 或 用 其 它 的 方法 


来 提高 码 环 的 鉴 相 精 度 。 国 外 学 者 在 快速 传 里 叶 变 换 的 基础 上 提出 内 插 的 思想 ， 提 取 相 


关 峰 中 的 最 高 点 与 最 低 点 ， 进 行 插值 处 理 提 取 相 位 ， 其 解决 了 侍 里 叶 变 换 计算 量 大 的 问题 ， 


但 是 精度 没有 提高 ; 国内 学 者 在 其 基础 上 提出 二 次 插值 的 方法 来 获取 更 精准 的 码 相位 测量 结 


果 ""， 结 果 表 明 在 信号 强 的 情况 下 有 良好 的 表现 ， 另 外 一 些 学 者 采用 最 小 二 乘法 拟 合 的 技 
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术 提高 码 相位 测量 精度 ， 该 方法 相 较 于 插值 法 有 实现 简单 、 精 度 高 的 特点 “， 但 其 测量 范围 


较 小 , 仅 适 合 码 相 位 误差 在 一 个 采样 点 以 内 测量 。 在 此 基础 上 有 学 者 提出 基于 整体 最 小 二 乘 


法 的 估计 算法 “， 虽 然 精度 有 所 提高 ， 但 需要 进行 矩阵 的 奇异 值 分 解 等 运算 ， 运 算 量 过 大 ; 


深 空 探测 中 使 用 再 生 伪 码 来 提高 码 相位 测量 精度  ， 再 生 伪 码 信号 的 码 相位 估计 依赖 码 元 


中 存在 显著 的 时 钟 分 量 , 利用 这 一 特性 ,依次 实现 码 片 时 钟 同步 和 码 相 位 同步 ， 两 者 实质 上 


分 别 对 应 小 数位 和 整数 位 的 码 相 位 值 ， 但 目前 没有 应 用 在 GNSS 接收 机 中 ; GNSS 接收 机 中 提 


二 | 


高 码 跟踪 环 相位 的 精度 方式 主要 还 是 在 环 路 的 参数 设置 与 选取 上 ,通过 设置 自 适 应 变化 


环 路 带宽 ， 实 时 调整 环 路 的 参数 , 来 提高 跟踪 环 路 精度 。 还 有 研究 提出 使 用 卡尔 曼 滤 波 跟踪 


环 代替 传统 码 相位 跟踪 环 ““， 其 提高 了 环 路 的 跟踪 性 能 与 动态 范围 ， 降 低 了 码 相 位 误差 与 


拌 动 , 但 是 计算 量 过 大 不 易于 实际 应 用 ; 文献 [16] 提 出 了 使 用 相位 条 纹 对 码 相 位 进行 测量 的 


方法 , 结果 表明 相位 条 纹 法 对 比 与 传统 GPS 接收 机 中 鉴 相 精度 有 所 提高 , 但 是 由 于 是 使 用 仿 


真 数据 进行 模拟 ， 缺 乏 实际 测量 数据 的 验证 。 


针对 以 上 问题 , 本 文 在 文献 [16] 的 基础 上 引入 实际 的 GPS 信号 , 针对 伪 码 信号 的 互相 关 功 


率 谱 ， 对 传统 的 跟踪 环 进行 改进 ， 提 出 一 种 高 精度 的 码 相 位 测量 方法 。 


1 相位 条 纹 法 测量 码 相 位 
1. 1 相位 条 纹 法 基本 原理 


相位 条 纹 方法 通常 用 于 测量 两 台 接收 机 之 间 的 距离 , 同一 信号 经 不 同 路 径 到 达 两 台 接 收 机 


后 , 通过 相关 计算 得 到 相位 差 , 从 而 计算 出 到 达 时 间 差 r, 乘 以 光速 得 到 两 台 接收 机 的 距离 。 


设 两 个 接收 机 接收 到 的 信号 分 别 为 (7) 和 (1)， 其 中 成 (1)= (1 一 +) ，7T 为 信号 到 达 


两 台 接 收 机 之 间 的 时 延 。 对 两 路 时 域 的 信号 做 传 里 叶 变 换 ， 并 进行 互相 关 计 算 , 得 到 的 互相 


关 功 率 谱 为 : 


Si,(®) = FP(0):F, (0)=F (Wer (1) 
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互相 关 功 率 谱 的 相位 : 
凶 Orz=27zr: 太 (2) 
可 以 看 出 , 在 两 个 接收 机 的 时 延 z 一 定 的 情况 下 ,， 互 功率 谱 相 位 乡 与 信号 的 频率 了 为 线性 


关系 ， 和 斜率 为 27z 。 通 常 相位 $ 的 取 值 范围 为 [-zx,7z]， 将 G~f 斜 线 按照 [-x, 如 的 范围 切 


割 成 周期 性 的 条 纹 ， 条 纹 之 间 的 间隔 为 7 ， 其 中 : 
二 


二 = (3) 
2NT 和 克 
而 条 纹 出 现 的 频率 三 : 
1 
二 (4) 
fo 
即 通过 统计 条 纹 出 现 的 频率 ， 可 得 到 信和 号 的 时 延 。 
4 才 1 1 1 1 | 1 1 1 1 | 1 h 1 1 1 1 1 1 | 1 h 1 1 由 


相位 由/rad 


T T T T T T T T T je 


0 1 2 频率 f/Hz 3 4 5 


图 1.， 相位 条 纹 示例 


Fig.1 Phase strip example 


1.2 相位 条 纹 法 测量 接收 机 码 相位 


在 GNSS 接收 机 当中 ， 接 收 的 信号 经 过 下 变频 剥离 载波 后 得 到 基带 信号 ， 基 带 信 号 与 本 地 


NC0 产生 的 伪 码 信号 可 以 视 作 是 两 路 相关 的 具有 时 延 的 信号 。 通 过 计算 出 其 中 的 待 测 时 延 可 


以 得 到 本 地 信号 与 接收 到 的 信号 的 码 相位 差 . 从 而 得 到 更 高 精度 的 码 相位 测量 值 。 接收 机 在 
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下 变频 后 通常 会 得 到 I、Q 两 个 正 交 的 支 路 ， 选 取 工 支 路 的 即时 信号 剥离 载波 并 与 本 地 NCO 


产生 的 码 信号 做 相关 运算 ， 得 到 互相 关 功 率 谱 ， 通 过 统计 出 现 的 相位 条 纹 频率 ,， 测 得 接收 信 


号 的 时 延 。 


但 由 于 接收 到 的 导航 信号 的 信 噪 比较 低 , 直线 拟 合 的 效果 较 差 , 故 首 先 对 接收 的 信号 与 本 


地 的 信号 进行 时 域 上 的 分 段 夺 加 ， 提 高 信号 的 信 吕 比 。 对 相位 条 纹 曲 线 G( 卫 ) 再 做 一 次 全 里 


叶 变 换 ， 当 条 纹 频 率 为 f。 ， 条 纹 间 隔 了 =1/ f， 单 个 条 纹 B (了) 可 以 表示 为 : 


1 1 
hn-1xlu1+ 去 | 中 去] (5) 


对 下 (用) 做 传 里 叶 变 换 ， 其 频谱 中 的 幅度 为 : 


sin A COS 


|@,Gx|=2x 一 、 (6) 


因为 条 纹 以 了 为 周期 出 现 ， 周 期 性 条 纹 的 相位 传 里 叶 变化 为 ; 


6 - 2 


Pw)|=2x > (7) 


n=]1 


当 式 (8) 中 的 n=1 


由 上 式 可 知 ， 相 关 条 纹 g( 了 ) 在 进行 傅 里 叶 变 化 后 得 到 的 幅度 谱 


时 ，|(w)| 可 得 到 最 大 值 ， 此 时 y=2x fh。 


男 外 ,由 于 接收 机 在 跟踪 信号 时 , 本 地 NC0 产生 的 伪 码 与 接收 到 的 伪 码 信号 小 于 1 个 码 


片 ， 时 延 z 较 小 ， 相 位 条 纹 的 出 现 频率 较 低 ， 为 了 更 加 明显 的 获得 相位 条 纹 ， 可 以 在 本 地 接 


收 信号 时 预先 补偿 一 个 时 延 r ， 令 m >>T， 在 整体 时 延 变 大 以 后 可 以 获得 频率 条 纹 数 量 ， 


有 利于 计算 出 更 为 精确 的 时 延 。 在 通过 条 纹 频 率 计算 出 z +z 以 后 ， 减 去 补偿 时 延 r ， 最 后 


得 到 待 测 时 延 z ， 即 可 得 到 码 相 位 测量 值 。 


相位 条 纹 法 测量 码 相 位 流程 如 图 所 示 : 
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$(7) 


Ya a 


本 地 复制 载波 


去 除 补 偿 时 延 


计算 互相 关 相 傅 里 叶 变 化 提 
位 谱 取 条 纹 频 率 


本 地 复制 码 


补偿 时 延 


NAC 


图 2， 相 位 条 纹 图 步 又 


Fig.2 Phase strip method steps 


1. 3 相位 条 纹 方法 提高 鉴 相 精 度 理论 分 析 


克拉 美 罗 下 界 CRLB (Cramer-Rao Lower Bound) 用 于 计算 无 偏 估计 中 能 够 获得 的 最 佳 估 


计 精 度 ， 因 此 经 常用 于 计算 理论 能 达到 的 最 佳 估 计 精 度 ， 和 评估 参数 估计 方法 的 性 能 。 下 面 


通过 分 析 相 位 条 弘法 中 的 CRLB。 
信号 经 过 下 变频 后 的 中 频 信号 为 : 


Xx(n)= A:D(t)cos[2z fyntT]+n, (8) 


其 中 A 为 信号 的 幅 值 ，7 为 待 测 码 相位 ，nj 为 方差 为 o 高 斯 白 噪声 ， 其 余 变量 均 为 常量 。 


其 似 然 函 数 为 : 


| x[n]— Acos(2z fin+7)] 


XT) = exp| 一 上 9 (9) 
BD 2x0” D 20°” 


求 其 对 数 一 阶 偏 导 ， 为 : 


4 N-l 
Oln pCG;7) = [xlnlsin(2z fin +7)—Asin(4x fon+ Se] (10) 


0 困 0 地 


其 二 阶 偏 导 为 : 


2 N-l 
- oe 9- > [xtnlcos(Qz fnt+7)—Acos(4zfon+27)] UL 
O 元 = 人 


二 阶 偏 时 的 期 望 为 : 


2 
如 In px;7) 
60 


其 中 期 望 E{X[A]) = Acos(2zfi; +7)， 上 式 期 望 化 简 得 : 


| 一 a SE {x[n]} cos(2z fin+7)— Acos(4z font+ 27) | (12) 
oO n=0 


OInpxn)| A 加 NA? 
-EF ;= 1— > cos(4x fon+27) | =NxSNR (13) 
| og 20” 2 之 oR 20° 


其 中 NN 为 信号 长 度 ，SNR 为 信号 的 信 噪 比 。 所 以 ， 中 频 信号 的 克 劳 美 罗 届 CRLB 为 : 


(14) 


Var{z} => 
NxSNR 


相位 条 纹 方 法 在 计算 时 通过 时 域 仅 加 的 方法 提高 了 信号 的 信 骂 比 , 因此 可 以 使 克拉 美 罗 界 更 


低 ， 提 高 估计 精度 。 


2 实际 GNSS 信号 实验 


本 次 实验 采用 数字 采集 卡 采集 的 GPS 信和 号 来 对 上 述 方法 进行 验证 。 选 用 的 GPS 采集 卡 为 


NSL Stereo， 其 采样 频率 为 26M， 中 频 载 波 频率 为 6.5MHz。 对 GPS 的 PRN32 卫星 数据 进行 


采集 ， 其 中 信号 长 度 为 100us， 码 速率 为 1.023MHz， 设 定 的 补偿 时 延 r 为 600ns。 后 续 的 处 


理 过 程 在 Matlab 上 完成 。 中 间 的 处 理 结 果 如 图 5~ 图 7 所 示 : 


| 


0.4 0.6 0.7 0.8 0.9 | 


0.5 
时 间 t/s x10% 
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图 3. 采集 的 中 频 信 号 


Fig3. IF signal collected 


相位 /rad 


图 4. 本 地 码 与 接收 到 的 码 互相 关 功 率 谱 相 位 条 纹 p( 了) 


Fig4. Local code and received code cross-correlation power spectrum phase strip of) 


0.8 


频率 wI0™® 


图 5. 对 相位 条 纹 进行 傅 里 叶 变化 提取 条 纹 频 率 (W ) 


Fig.5 Perform Fourier transformation on phase strip to extract Strip frequency DY) 


5 为 采集 的 数据 中 截取 的 100us 的 长 度 部 分 , 通过 对 其 混 频 处 理 后 剥离 中 频 载波 , 留 下 


chinaXiv:202107.00042vVv1 


ChinaXiv 合 作 期 刊 


伪 码 信号 ， 与 本 地 产生 的 PRN32 号 卫星 的 码 做 相关 运算 ， 得 到 其 相关 频率 条 纹 G(f)， 如 图 


6 所 示 。 对 G(f) 再 进行 伟 里 叶 变 化 ， 其 得 到 的 功率 谱 的 |®(w)| 的 最 大 值 对 应 的 频率 即 为 所 


需要 测 得 的 时 延 z+zo ， 即 图 7 中 |@®(w)| 最 大 值 对 应 的 频率 为 所 测 时 延 t+T6。 经 计算 ， 


5 十 70=607.6ns, 减 去 ro =600ns, 最 后 得 到 7 为 =7.6ns。 乘 以 光速 后 得 到 的 鉴 相 误差 为 2.28m。 


接收 机 的 跟踪 环 路 在 进行 稳定 的 跟踪 后 ， 其 理论 热 噪声 误差 的 计算 如 下 : 


) (15) 


将 误差 转化 为 距离 ， 可 得 : 


4 | 也 1 
Om = 1+ (16) 
”2x\CNV, 2CNT 


在 公式 (9) 中 ，CN。 为 接收 机 的 载 品 比 ，B， 为 噪声 带宽 ,T 为 相关 积分 时 间 ， 在 本 次 实 


际 采 样 的 数据 中 ， 载 噪 比 通过 采样 数据 计算 出 为 CN。 =42. 6dB-Hz。 传 统 的 接收 机 积分 时 间 


为 T=lms，B, =2Hz， 在 载 品 比 为 42. 6dB-Hz 时 ， 其 理论 热 噪 声 误差 上 限 为 0 w =36. 059m， 


实际 使 用 传统 延 时 锁定 环 测量 结果 为 5. 39m。 通 过 相位 条 纹 计算 出 的 误差 为 2. 28m， 优 于 传 


统 锁 相 环 中 非 相 干 点 积 功率 法 鉴 相 的 理论 误差 。 


4. 结 论 


传统 GNSS 接收 机 在 实际 的 跟踪 中 的 鉴 相 误差 在 0. 1 个 码 片 左右 ， 即 为 2. 93m， 通 过 本 文 


中 的 方法 得 到 的 误差 为 2. 28m， 人 性 能 提高 了 23%， 在 实际 的 接收 机 当中 ， 为 提供 高 精度 的 码 


相位 计算 给 出 了 一 种 方法 ， 从 而 得 到 了 更 高 精度 的 伪 距 测量 值 。 而 在 后 续 的 工作 中 ， 需 要 对 


相位 条 纹 方法 中 的 参数 的 选取 进行 进一步 的 研究 。 


A high-preclsion code phase measurement method based 


on phase strip 


Wenjie Sun!“，Zhaorui Wang' 
(1. Chinese Academy of Sciences, National Astronomical Observatories of China, Beijing 100012，China 
2. ”University of Chinese Academy of Sciences，School of Astronomy and Space Science, Beijing 100049, China) 

Abstract: After the GNSS (Global Navigation Satellite System) receiver receives the satellite 
signal, it acquires and tracks the two basic measurement values of code phase and carrier phase. 
Calculate the distance between the satellite and the receiver based on the measured value. The 
carrier phase measurement value has a higher measurement accuracy than the code phase 
measurement value, but because the carrier phase measurement value has the problem of 
calculating the ambiguity of the whole cycle, in the actual, the high-precision phase measurement 
result requires more calculation cost. Therefore, in actual calculations, obtaining a 
higher-precision code phase measurement result has Practical measurement significance. Based on 
the phase strip technology, this paper proposes a higher-precision code phase measurement 
method by extracting the cross-correlation power Spectrum to obtain a higher-precision code phase 
measurement value, thereby obtaining a higher-precision code pseudo range measurement result, 
And verify the effectiveness of this method through actual GNSS signals. 
Key Words: GNSS; Phase Strip; Code Phase 
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